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One of the key issues in cosmology is to establish the nature of dark energy, and to determine 
whether the equation of state evolves with time. When estimating this from distance measurements 
there is a degeneracy with the matter density. We show that there exists a simple function of the 
dark energy equation of state and its first derivative which is independent of this degeneracy at all 
redshifts, and so is a much more robust determinant of the evolution of dark energy than just its 
derivative. We show that this function can be well determined at low redshift from supernovae using 
Gaussian Processes, and that this method is far superior to a variety of parameterisations which are 
also subject to priors on the matter density. This shows that parametrised models give very biased 
constraints on the evolution of dark energy. 


I. INTRODUCTION. 


Since the discovery in 1998 that the Universe is under¬ 
going a period of accelerated expansion 0,i , consider¬ 
able work has been done trying to understand what it 
is driving the phenomenon. The standard explanation 
beyond a simple cosmological constant relies on intro¬ 
ducing a new cosmic fluid, dubbed dark energy, which 
possesses negative pressure in order to accelerate the ex¬ 
pansion. Dark energy is parametrised via an equation of 
state w{z) = where pde is its pressure, pde is its 
energy density and w{z) stands for a possible dependence 
of the equation of state with the redshift z. The standard 
case is given by w(z) = —1, which represents the cosmo¬ 
logical constant A, where together with cold dark matter 
(CDM) form the cosmic concordance model ACDM. 

Extensive observational efforts have been carried out 
in order to falsify w(z) = —1 through several observables, 
as type la supernovae (SNe la), temperature anisotropies 
of the cosmic microwave background (CMB) and baryon 
acoustic oscillations (BAOs). As there is no compelling 
model to explain the acceleration, it is often simplest to 
use simple parameterisations for w, where a dependence 
with the redshift can be accommodated by adding a new 
parameter. On the other hand, an interesting and inde¬ 
pendent way to test for deviations of the ACDM model is 
to use non-parametric methods, which are suitable for de¬ 
scribing how w(z) evolves not restricted by simple para¬ 
metric models. Examples considered in the literature in¬ 
clude the principal component analysis method (PCAs) 
[1], gaussian processes (GPs) d H and local regression 
smoothing Q. 

Reconstructions have been limited by degeneracies in- 
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volving the matter density parameter and the cur¬ 
vature density parameter From one side, flk can 
be inferred due its geometrical structure [ 3 , whereas for 
flm there is a so-called dark degeneracy d whereby dark 
energy and dark matter are hard to separate. 

Here we show that there is a simple function of w(z) 
and its first derivative which is completely independent 
of flm- This function can be used to look for a time 
evolution of the dark energy equation of state without 
any dark matter degeneracy. We use data from SNIa 
and i?(z) measurements from cosmic chronometers and 
BAOs to constrain this new function, along with the evo¬ 
lution of w. We compare parametric reconstructions to 
a GP reconstruction and show that GPs provide much 
stronger constraints than parameterisations at low red- 
shifts and are broadly consistent with AGDM. 


II. THE EVOLUTION OF DARK ENERGY AND 
THE DISTANCE 


The dimensionless comoving distance D(z) can be 
written as a function of the luminosity distance dL{z) 
as 


D{z) = {Hq/c){ 1 + z) '^dL{z), (1) 


where Hq is the Hubble constant. The effective dark 
energy equation of state is 


2{1 + z)D" + 3D 
30^(1-bz)3i:)'3-3D'’ 


where a prime denotes differentiation with respect to the 
redshift. Throughout the paper we assume Utk — 0, al¬ 
though we note a degeneracy between w and flk can be 
important at high redshifts Given distance-redshift 
data, we can determine w from this relation. However, 
the value of Dm must be provided a priori from indepen¬ 
dent methods in order to test for deviations of lu = — 1 
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M- Similarly, one can write the effective adiabatic 
sound speed for dark energy as 


1 ?// 

2 (1 + z)D'D"' - 2D'D" - 3(1 + z)D”^ 

3 D' [2D" + 312^(1 + zfD'^] ’ 


Combining with w{z) we define 


F{z) = 


l + w(z) 2 


w{z) 




( 5 ) 


= + + ^ + w{z) 

3 w{z) 


2 (1 + z) [(1 + z)D'D"' - 2D'D" - 3(1 + z)D"‘^] 

3 D' [2{l + z)D'’ + ZD'] 


Note that F{z) is completely independent of Vim and is 
a function of the observable quantity D and its deriva¬ 
tives. Model-independent constraints on F - and con¬ 
sequently w - can be derived through non-parametric 
methods. Any deviation from F — 0 would rule out 
the flat ACDM model, and so forms a useful null test 
of concordance cosmology m (see [H for another null 
test independent of Vim)- Another way to think of the 
function F is to say that in the possible function space 
of w{z), F selects the subspace which is independent of 
Vim, while giving a measure of the first derivative of w. 


A. Constraining w and F 


B. Data 

For both the non-parametric and parametric analyses 
we use 580 luminosity distance measurements from SNe 
la 0 , where we included all the systematic errors, and 
26 H{z) measurements from cosmic chronometers [T^ 
and BAOs [l^. For the Hubble constant we adopt FIq = 
70.4±2.5kms-i MpC , which is in agreement with 
the revisited local value based on the NGC 4258 maser 
distance obtained by Efstathiou [0: FIq = 70.6±3.3 km 
s“^ Mpc“^ . 

As we assume spatial flatness {Vlk = 0), the H{z) mea¬ 
surements can be directly used to obtain the first deriva¬ 
tive of D: D'{z) = which allows us to incorporate 
the derivative of D as data in our GPs analysis. 

For the parametric analyses, we obtain the posterior 
probabilities (see e.g. [0) using emcee [l^, an affine- 
invariant ensemble sampler for Markov Chain Monte 
Carlo (MCMC) [2l|. For example, in the GPL case we 
use the posteriors of wq — Wa to derive constraints to 
w{z) sampling from these parameters, the same for F 
using the second line of Eq. [5l 

An important issue in our analysis is that distance 
moduli of SNe la were used already fitted for some nui¬ 
sance parameters in a given cosmological model [T3| . 
Therefore, a more model-independent procedure would 
be to introduce these nuisance parameters in our GP re¬ 
gression. On the other hand, such approximation was 
also taken in our parametric analyses, which makes a 
comparison between them consistent. 


We shall now constrain F in 2 different ways. Given 
SNe la and H{z) data, we need to find D{z) and its first 
3 derivatives. First we will use Gaussian Processes (GPs) 
which are a robust smoothing technique which allows for 
differentiation of data - this technique is summarized in 
the Appendix. We will use the third line in Eq. [S] to 
obtain constraints on F and Eq. [^with a suitable prior 
on Vim to derive constraints for w. Then we will compare 
these results with two parametric models. We choose 
a standard one and a new one which prodnces similar 
constraints on w but rather different constraints on the 
evolution. 

• A standard series expansion in 1 — a: w{z) = 

Wo ± generally called the GPL (Ghevalier- 

Polarski-Linder) parameterisation [0. 

• Gonstant F: wiz) = - ■ This 

model is constructed such that F(z) = 1 — C is 
constant. Note that (7=1 gives F = 0, which 
means that although F 0 rules out flat ACDM, 
F = 0 is also compatible with different models. 

Although the second of these may look rather unusual, 
these are actually rather similar in terms of their basic 
assumptions because the GPL parameterisation is con¬ 
structed so that w^a is constant. 


C. Results. 

In Fig. [T] we show constraints on w{z), and F{z) for 
the different methods. The parametrised models show 
broadly similar behaviour for the constraints on w{z) and 
give similar errors compared to those found using GPs, 
at least for low to moderate z (strong constraints at high 
z only come from the choice of parameterisation). In 
this case, the good constraints with GPs are obtained 
because of the prior adopted on Vim, where much weaker 
ones would be derived considering a broad flat prior on 
Vim - this arises simply from ([2]). For F(z) the situation 
is rather different, with very different constraints found 
for the different approaches. Constraints at all redshifts 
between parametrised models versus GPs are very differ¬ 
ent. With GPs F(z) is well constrained at small redshifts 
(z < 0.5) and weakly constrained for redshifts above that. 
Also it is interesting to point ont that for GPs the con¬ 
straints on F(z) are weakly dependent on the H{z) data, 
where they are slightly enlarged if we do not use this 
sample. In redshift regions well sampled by data good 
constraints are derived, while redshift regions sparsely 
measured are not able to distinguish between different 
values for F, when combined with the fact that changes 
in w affect w at low z most strongly. This is in stark 
contrast with the results from the parametric methods. 
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FIG. 1: w{z) constraints, top and corresponding F{z) constraints, bottom. The contours represent 68% and 95% confidence 
intervals. The left panels show the GP reconstruction. The middle panels depict the constraints the GPL parameterisation, 
while the right panels the constant F parameterisation. A gaussian prior on Ilm = 0.28 ± 0.02 is assumed for w, no prior is 
needed for F (though for the parametric models, the contours narrow slightly with a prior on Ilm). 








FIG. 2: Posterior for w(0), w'(0) and F{0) found using Gaussian Processes, versus parametric approaches. The top row 
has a fiat prior on Qrn, the bottom row has a Gaussian prior flm = 0.28 ± 0.02. For all the parametric approaches the 
constraints are strongly dependent on the prior. For GPs, a prior is necessary for w(0), w'(0) but not F. For parametric 
approaches, constraints on t(j(0) are reasonably consistent but achieve very different results for w'(0) and F’(O), effectively 
giving no meaningful constraints on the evolution of the EOS at all. GPs on the other hand give the strongest constraints on 
all variables when a prior on flrn is included (it was shown that the errors on w using this method are robust in [^1. Strong 
constraints on F are found using GPs independently of Qm- 
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which provide better results at high redshifts due to the 
imposed shape of the function a priori. 

In Fig. [5] we show the posteriors at z = 0 for all the 
different cases. With no prior on 17^ only F found by 
GPs has strong constraints, while individual constraints 
on w and w' are significantly weaker. With a prior on rim, 
the constraints on w{0), w'{0) and F{0) vary significantly. 
In particular, the parametric approaches favour both w 
and w' in different ways (CPL favours w nearer 0 and 
w' more negative, and the constant F model pulls in 
the opposite directions), while GPs give tight constraints 
approaching the union of the parameterised models. 

It is important to stress one needs to be careful of not 
deriving biased results due to a given value for ftm- For 
example, there is a tension for w'(0) when using GPs and 
our (somewhat arbitrary) choice of prior on Qm which 
completely disappears when using Planck values for 
rim and F[q: w'( 0) = —0.03 ± 0.23 (Icr). The H{z) data 
prefer lower values of iJo as it was shown in (2^. There¬ 
fore, GP reconstructed F gives strong constraints favour¬ 
ing a ACDM model, with F providing an unbiased, more 
powerful tool to look for deviations of the concordance 
model as a time-varying dark energy. 

D. Validating the method 

In order to analyse the robustness of our method, we 
performed 1000 simulations of SNe la data in a ACDM 
model with rim = 0.3, Fl^ = 73.8 km s“^ Mpc“^ and 
the same redshift distribution of the Union2.1 sample. 
For the sake of simplicity, we restricted ourselves only to 
SNe la data and we show the results for F{z) with an 
emphasis on F'(O). 

Fig 13] presents the average posterior for F given the 
1000 simulations, where the left panel has the result for 
F(0) and the right panel the result for F{z) for all red- 
shifts. As we can see, the width of the curve is consis¬ 
tent with the errors obtained, and there is a small bias 
much smaller than the errors. We also have checked that 
the results are not influenced by the fiducial cosmologi¬ 
cal model, since changing for XCDM or a time-varying 
equation of state for dark energy gave consistent results, 
as it is shown in Fig. 4. Bigger samples as the ones com¬ 
ing from the ongoing Dark Energy Survey will demand a 
thorough scrutiny of the bias, although it can be safely 
neglected for the current size of SNe la samples. 

III. CONCLUSIONS 

We have shown that for a flat model there exists a 
combination of w and its first derivative w' which can 
be found purely from Hubble normalised distance data 
independently of rim- This provides a route by which 
we can place constraints on the time evolution of w 
without requiring a prior on rim- At z = 0 this cor¬ 
responds to roughly the lines in the wq — Wa plane where 


Wa/3wo -I-1 -b Wo =const. We have shown that using GPs 
strong constraints can be placed on F{0) using current 
data, and with no assumptions on rim or w{z). This is to 
be seen in stark contrast to parametric models which have 
significantly weaker constraints, but an order of magni¬ 
tude. Furthermore, the constraints on F are similarly 
much stronger than constraints on either w or w' at the 
origin. Consequently, using E as a complimentary func¬ 
tion to w{z) might be a more useful way to characterize 
the time evolution of w than is currently done. 

Of course, if one wants w{z) explicitly the degeneracy 
with rim is unavoidable. In this sense, have we really 
dodged the degeneracy? We have presented a new per¬ 
spective on the question: since w{z) is usually a place¬ 
holder for whatever is ‘not-A’, it may well be a sub- 
optimal measure of ‘not-A’. A set of observables such 
as [F{z), F'(z), F''{z ),...], formed from ([5]) rephrase the 
A question as: do [F{z), F'{z), F"{z),.. .]z=o = 0? but 
now independently of rim- 
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Appendix A: Gaussian processes 

A Gaussian process is defined by a collection of ran¬ 
dom variables, any finite number of which have (consis¬ 
tent) joint gaussian distributions [2^. Basically speak¬ 
ing, while a gaussian variable is a distribution over ran¬ 
dom variables, a gaussian process is a distributions over 
random functions, characterized by a mean and a covari¬ 
ance function. 

In order to reconstruct a given function, we assume 
that this function is a realization of a GP with a given 
mean and covariance function. Assume we have N obser¬ 
vational data points y = {yi,..., y^) of a function f(z) 
at redshifts Z = {zi,..., zn). The errors of the observa¬ 
tions are given in the covariance matrix C. Additionally 
we have N observations of the first derivative f’{z), given 
by Z = (fi,..., 5 jy), y' = {y[, ■ ■ ■, y'ff) and C. Thus, the 
goal is to reconstruct the function which underlies the 
data at redshifts Z* = {zl,...,z%). We denote these 
function points as f* = (f(z ^),..., f{z’^)) and the nth 
derivative of the function at these redshifts as 

Generally, we incorporate the desired properties of the 
function under study in the covariance function. In this 
work we adopt the Matern(n = 9/2) covariance function 
(see [ 2 ^ for a discussion on the impact of the covariance 
function) 
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FIG. 3: Average posterior of one thousand simulations of SNe la data for a ACDM model with flm = 0.3 and Hq = 73.8 km 
s”*^ Mpc~^. The same redshift distribution of the Union2.1 sample was used. In the left panel is shown the behaviour of F 
at redshift zero, while in the right panel for the whole redshift range. The size of errors is fully compatible with the observed 
data, apart from a small bias which can be neglected for the current sample size, but that should be better characterized for 
bigger samples. 



FIG. 4: Typical realization of SNe la data with different fiducial models. The left panel uses a fiducial model with constant 
w = 1.1 which gives F = 0.1. The right panel uses a fiducial model with a time-varying equation of state with wq = 1 and 
Wa = 0.1. The solid line shows the value of the expected model. 
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where we have two hyperparameters aj and £ which de¬ 
scribe how the function changes in the y and x axis, re¬ 
spectively. 


The probability distribution for the points we want to reconstruct is gaussian with the mean and covariance given 
by 


/*(") = 




'K{Z,Z) + C 
K^^’°'>{Z,Z) 


K^^P{Z,Z) 

kP^\Z,Z)+C_ 
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and 


cov 



Z*)-{K<-^’°\Z*, Z) {Z*, Z)) 


'K{Z,Z) + C X(0’1)(Z,Z) 
K^^'°\Z,Z) K^^-'^\Z, Z) + C_ 




(A3) 

Here, the covariances are written as matrices K{X,X) with {K{X, X)}ij = k{xi,Xj). denotes the nth 

derivative of k with respect to the first argument and the mth derivative with respect to the second argument. 

The hyperparameters are determined by maximizing the likelihood 


ln£ 


\np{y\Z,Z,af,e) 

K{Z,Z) + C K^°'^\Z,Z) 
K^^’^\Z,Z) K^^’^\Z,Z) + C 


{y^ y'l 


K{Z,Z) + C X(0’b(Z,Z) 



-bl- 

V»7 

2 


(A4) 


7V + 

2 


In 27r. 


While the best approach would be to marginalize over the hyperparameters, it was shown in that for a sample 
of the same size considered in this work, maximization provides basically indistinguishable results at low redshifts, 
which is the region where good constraints were derived. 

We use GaPP (Gaussian Processes in Python) to derive the GP results. The code is free to use and was applied 
in many situations [^ . 
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